Phase transition in the 3-state Potts antiferromagnet on the diced lattice 
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We prove that the 3-state Potts antiferromagnet on the diced lattice (dual of the Kagome lattice) 
has entropically-driven long-range order at low temperatures (including zero). We then present 
Monte Carlo simulations, using a cluster algorithm, of the 3-state and 4-state models. The 3- 
state model has a phase transition to the high-temperature disordered phase at v = e J — 1 = 
—0.860599 ± 0.000004 that appears to be in the universality class of the 3-state Potts ferromagnet. 
The 4-state model is disordered throughout the physical region, including at zero temperature. 
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The g-state Potts model 0, Q plays an important 
role in the theory of critical phenomena, especially in 
two dimensions [!, 0, @| , and has applications to various 
condensed-matter systems . Ferromagnetic Potts mod- 
els are by now fairly well understood, thanks to universal- 
ity; but the behavior of antiferromagnetic Potts models 
depends strongly on the microscopic lattice structure, so 
that many basic questions must be investigated case-by- 
case: Is there a phase transition at finite temperature, 
and if so, of what order? What is the nature of the low- 
temperature phase(s)? If there is a critical point, what 
are the critical exponents and the universality classes? 
Can these exponents be understood (for two-dimensional 
models) in terms of conformal field theory Q? 

One expects that for each lattice C there exists a value 
q c {£) such that for q > q c {C) the model has exponential 
decay of correlations uniformly at all temperatures, in- 
cluding zero temperature, while for q = q c (C) the model 
has a zero-temperature critical point. For q < q c (C) any 
behavior is possible; often (though not always) the model 
has a phase transition at nonzero temperature, which 
may be of either first or second order Q. The first task, 
for any lattice, is thus to determine q c . 

Some two-dimensional antiferromagnetic models at 
zero temperature have the remarkable property that 
they can be mapped onto a "height" (or "interface" or 
"SOS-type") model Experience tells us that when 
such a representation exists, the corresponding zero- 
temperature spin model is most often critical [9j|. The 
long-distance behavior is then that of a massless Gaus- 
sian with some (a priori unknown) "stiffness" K > 0. 
The critical operators can be identified via the height 
mapping, and the corresponding critical exponents can 
be predicted in terms of the single parameter K . Height 
representations thus provide a means for recovering a sort 
of universality for some (but not all) antiferromagnetic 
models and for understanding their critical behavior in 
terms of conformal field theory [|| . 



In particular, when the g-state zero-temperature Potts 
antiferromagnet on a two-dimensional lattice C admits a 
height representation, one expects that q = q c {C)- This 
prediction is confirmed in all heretofore-studied cases: 3- 
state square-lattice lil. [Tol 11, 12 1, 3-state Kagome 13, 



14| , 4-state triangular [15[ , and 4-state on the line graph 
(= covering lattice) of the square lattice 14, 1( 



We now wish to observe that the height mapping 
employed for the 3-state Potts antiferromagnet on the 
square lattice Q carries over unchanged to any planar 
lattice in which all the faces are quadrilaterals. One 
therefore expects that q c = 3 for every (periodic) plane 
quadrangulation. 

The diced lattice (Fig. [1]) is a periodic tiling of the 
plane by rhombi having 60° and 120° interior angles; in 
particular, it is a plane quadrangulation in which all ver- 
tices have degree 3 or 6. The diced lattice is the dual of 
the Kagome lattice, which is in turn the medial graph of 
the triangular and hexagonal lattices. 

In this Letter we give a rigorous proof that the 3-state 
diced-lattice Potts antiferromagnet has a phase transi- 
tion at nonzero temperature. This shows that, contrary 
to theoretical expectations, g c (diced) > 3. It provides, 
moreover, the first example of a bipartite two-dimensional 
lattice in which q c > 3 (but see below). 

Furthermore, this model provides a simple concrete ex- 
ample of entropically-driven long-range order, in which 
coexistence between regions of different types of order on 
one sublattice is disfavored because it reduces the free- 
dom of choice of spins on the other sublattice. Though 
this idea is intuitively appealing, it is difficult to deter- 
mine, in any specific case, whether the mechanism is 
strong enough to produce long-range order. Here we are 
able to resolve this question rigorously by a conceptually 
simple though numerically delicate Peierls argument. 

The existence (though not the nature) of a phase 
transition in this model is not, however, totally unex- 
pected. A decade ago, Jensen et al. jTz} computed 
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Figure 1: A diced lattice of size 3x3 with periodic boundary 
conditions (edges depicted with thick black lines). The full 
circles show the sites of degree 6, which form a triangular 
lattice (edges depicted with thin gray lines). The open circles 
show the sites of degree 3, which form a hexagonal lattice 
(edges depicted with thin dashed green lines) that is the dual 
of the triangular lattice. Periodic boundary conditions are 
implemented by identifying border sites with the same label. 



low-tempcraturc expansions for the 3-state and 4-state 
Potts models on the Kagome lattice and found, among 
other things, indications of singularities in the unphysi- 
cal region at v = -3.486 ± 0.003 and v = -3.38 ± 0.06, 
respectively (here v = e J — 1 where J is the nearest- 
neighbor coupling and we take (3 = 1). Shortly there- 
after, Feldmann et al. [l8| used the duality v <— > q/v 
of g-state Potts models on planar lattices to deduce pre- 
dictions for the singularities of the 3-state and 4-state 
Potts models on the diced lattice: v = -0.8607 ± 0.0008 
and v = —1.18 ± 0.02, respectively. The latter occurs 
in the unphysical region at v < —1. suggesting that the 
4-state diced-lattice antiferromagnet lies in a disordered 
phase at all temperatures, including zero temperature. 
The former, by contrast, lies within the physical anti- 
ferromagnetic regime at J = —1.971 ± 0.006. If these 
predictions are correct, we have 3 < g c (diced) < 4; crude 
linear interpolation suggests q c (diced) ks 3.4. 

One might, however, worry that the errors in the se- 
ries extrapolation are radically larger than estimated and 
that the diced-lattice singularity lies not at v f» —0.86 
but instead at the theoretically expected v = —1 (which 
is not, after all, so far away). It is thus important to 
obtain independent evidence on the location of the phase 
transition (if any) in the 3-state diced-lattice Potts an- 
tiferromagnet. We do this in two steps: a mathemati- 
cally rigorous proof of a phase transition at nonzero tem- 
perature, which furthermore elucidates the entropically- 
driven nature of the corresponding long-range order; and 
Monte Carlo simulations to locate the phase transition 
and investigate its properties. 



Proof of phase transition. We shall prove that, at all 
sufficiently low temperatures, there is antifcrromagnctic 
long-range order in which the spins on the triangular 
sublattice take preferentially one value and the spins on 
the hexagonal sublattice take more-or-less randomly the 
other two values. To implement the heuristic idea of 
entropically-driven order, we evaluate the cost of coexis- 
tence between regions of unequal spins on the triangular 
sublattice by integrating out the spins on the hexagonal 
sublattice, yielding a g-state Potts model on the trian- 
gular lattice with 3-body interactions on the triangles, 
namely, Boltzmann weights (wi, W2, W3) = (q+3v+3v 2 + 
v 3 , q + 3v + v 2 , q + 3v) according as the triangle has 1, 2 or 
3 distinct spin values [l|| . We then apply a Peierls argu- 
ment to this general triangular-lattice model and prove 
that there exists ferromagnetic long-range order in an 
open region of (w^/wi, u^/w^-space that includes the 
point (w2/wi,W3/wi) = (1/2,0) corresponding to q = 3 
at zero temperature (v = — 1). 

Choose a large box A and fix all boundary spins in the 
same state. Given a spin configuration on the triangular 
lattice, draw Peierls contours on the dual hexagonal lat- 
tice in the usual way, i.e., separating unequal spin values 
on the triangular lattice. We thus obtain a spanning sub- 
graph of the hexagonal lattice in which the vertices have 
degree 0, 2 or 3 according as the corresponding triangle 
has 1, 2 or 3 distinct spin values. 

Consider first the case W3 = (this covers the zero- 
temperature 3-state diced-lattice model, i.e., q = 3 and 
v = —1). Then the Peierls contours have no vertices 
of degree 3, so they are disjoint unions of self-avoiding 
polygons (SAPs) on the hexagonal lattice. Each contour 
edge gets a weight W2 jw\ , and each contour gets an addi- 
tional weight q~ 1 to count the possible values for the spin 
change modulo q when crossing the contour. If the prob- 
ability of having at least one contour is less than (q—l)/q, 
then the spin at the origin has probability greater than 
1/q of being in the same state as the boundary condition, 
hence there is long-range order. This occurs whenever 
00 

(tito/iai)" < 1/q, (1) 

n— 6 

where q^ is the number of n-step hexagonal-lattice 
SAPs surrounding the origin of the triangular lattice, or 
equivalently the first area-weighted moment for n-step 
hexagonal-lattice SAPs modulo translation. To bound 
this sum, we use the exact values of g„ for 6 < n < 140 
[U and the bound q£ ] < (n 2 /36) 1.868832"- 2 for even 
n > 142 [2l]]. For q = 3 we deduce long-range order 
whenever W2/W1 < 0.503417, which barely includes the 
desired value W2/W1 = 1/2 [22|. 

The case W3 > is more complicated because of the 
presence of degree-3 vertices, but it can be shown [23| 
that if the Peierls inequality holds at a given value of 
W2/W1 when W3 = 0, then it will also hold at that same 
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Figure 2: Coarse plot for the Binder ratio 7?. Dotted vertical 
line marks the critical point predicted in fvH , Il8| ] . Curves are 
straight lines connecting points, meant only to guide the eye. 



value of W2/W1 whenever w^/wi is sufficiently small. This 
proves that the 3-state diced-lattice antiferromagnet has 
long-range order on the triangular sublattice at all suffi- 
ciently low temperatures. With a little more work [23j |. 
we can prove antifcrromagnctic long-range order on the 
whole diced lattice. 

Monte Carlo simulation. We simulated the diced- 
lattice Potts antifcrromagncts for q = 2,3,4 on L x L 
lattices (3 < L < 768) with periodic boundary condi- 
tions, using the Wang-Swcndsen-Kotecky (WSK) clus- 



ter algorithm [24(. Since the diced lattice is bipartite, the 
WSK algorithm is guaranteed to be ergodic [3, [25| and 
there is reason to hope that critical slowing-down might 
be absent (as for the square lattice [ID]) or at least small. 

We measured the energy £, the sublattice magnetiza- 
tions jVJhex and Mtri-, and the second-moment correla- 
tion length £ [26[ . We focussed attention on the Binder- 
type ratio R = (A^ tagg ) / '(-^tagg) 2 where -M st a gg = 
■Mtri — .Mhex, which tends in the infinite-volume limit to 
/(<?— 1) in a disordered phase and to 1 in an ordered 
phase, and is therefore diagnostic of a phase transition. 
The ratio £/£ plays a similar role. Finally, we studied 
(A^stagg) in order to estimate the leading magnetic crit- 
ical exponent. 

The Ising (q = 2) data confirm the exact solution 
For q — 4, we find a finite correlation length uniformly 
down to zero temperature, with "f~ 1-85 as v I — 1. 

For q = 3, a plot of R versus v for lattice sizes 3 < L < 
96 is shown in Fig. [2l and shows a clear order-disorder 
transition at v c w —0.86. Finer plots of R and £/L near 
the transition, for 48 < L < 768, are shown in Figs. [3] 
and [4] We fit the data to Ansiitze obtained from 
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Figure 3: Fine plot for the Binder ratio R. Curves are our 
fits to ([2]). Symbol ^ indicates estimates of v c and R c . 
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Figure 4: Fine plot for £/L. Curves are our fits to 
Symbol ^ indicates estimates of v c and (£/L) c . 
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O c + ai(v - v c )L^ u + a 2 (v - v c ) 2 L 2 ^ 



by omitting various subsets of terms, and we systemati- 
cally varied L m in (the smallest L value included in the 
fit). We also made analogous fits for (A4g tagg ) jU 1 ^ . 
Comparing all these fits, we estimate the critical point 
V c = -0.860599 ± 0.000004, the critical exponents v = 
0.81 ± 0.02 and j/v = 1.737 ± 0.004, and the univer- 
sal amplitude ratios R c = 1.170 ± 0.007 and (£/L) c = 
0.995 ± 0.007 (68% subjective confidence intervals, in- 
cluding both statistical error and estimated systematic 
error due to unincluded corrections to scaling). These 
exponents are in excellent agreement with the values 
for the 3-state Potts ferromagnet: v = 5/6 w 0.833, 
7/V = 26/15 w 1.733 and R c = 1.1711±0.0005 0. This 
confirms our expectation that the 3-body-intcracting 
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triangular- lattice ferromagnet, obtained by integrating 
out the hexagonal sublattice, lies in the universality class 
of the 3-state Potts ferromagnet. 

For the triangular-sublatticc spontaneous magnetiza- 
tion M , defined by M$ = lim i _K X) (Al^ i )/^ i , we find 
M = 0.936395 ± 0.000006 at v = -1, which is not 
far from the Peierls bound M > 0.90497 [HI and 
the heuristic estimates Mo ~ 125/128 ps 0.97656 and 
M 21/22 ps 0.95455 

Details of the simulations will be reported later [23| ■ 

Final remarks. Here is another construction that pro- 
duces bipartite planar lattices (not, however, plane quad- 
rangulations) with q c > 3 and indeed with q c arbitrarily 
large. Let £ be any lattice, and let £2 be the lattice ob- 
tained from £ by subdividing each edge into two edges 
in series. Then the Potts series law (ui, V2) >— ► ^1^2/(9 + 
vi + V2) [IH implies that the Potts model on £2 has a 
phase transition whenever v 2 /(q + 2v) = i> cr it, ferro, c(q)- 
In particular, if u C rit,ferro,£(?) < 1/(9 — 2), then there 
is a solution v e (—1,0), so that q c {£2) > 9- For in- 
stance, the triangular lattice T has a ferromagnetic criti- 
cal point when v 3 + 3u 2 — 9 = 0, from which we conclude 
that 9 C (T 2 ) ps 3.117689. Furthermore, Wicrman [3(j has 
constructed periodic plane triangulations T(k) [obtained 
by subdividing triangles in the triangular lattice] whose 
bond percolation thresholds tend to zero as k — ► 00; and 
the Potts series-parallel laws show, more generally, that 
for each 9 > 1 one has lim^^ v cri t,ferro,T(fc) (9) = [31 1- 
It follows that linifc^oo q c (T{k)2) = +00. 
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